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We explore the influence of external perturbations on the energy levels of a Hamiltonian drawn 
at random from the Gaussian unitary distribution of Hermitian matrices. By deriving the joint 
distribution function of eigenvalues, we obtain the (n,m)-point parametric correlation function of 
the initial and final density of states for perturbations of arbitrary rank and strength. A further 
generalization of these results allows for the incorporation of short-range spatial correlations in 
diffusive as well as ballistic chaotic structures. 
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PACS numbers: 71.23.-k,05.45.Mt,73.21.-b 

Beginning with early pioneering work on resonances 
in complex atomic nuclei the statistical properties of 
general chaotic and disordered quantum structures have 
come under intense scrutiny ||^. In such systems, it 
is conjectured (and, in the case of disordered systems, 
proven explicitly [^) that the low-energy correlations of 
eigenvalues (and cigcnfunctions) of an individual system 
can be inferred from the ensemble average properties of a 
random matrix distribution whose members exhibit the 
same fundamental symmetries Q . For an N x N matrix 
Hamiltonian H drawn from, say, a Gaussian distribution, 
the statistical properties of the eigenvalues {e^} are com- 
pletely specified by the joint distribution function 



P({e.})oc|A(M)|'^e-(^/2)5:£i< 



(1) 



where A({ei}) — Yli^j (ci — fj), and /3 = 1, 2 and 4 for 
ensembles of orthogonal, unitary and symplectic symme- 
try, respectively. From this result, it follows that the 
n-point correlator of density of states (DoS) Rn{{^a}) = 
(Y[a=i ^ (^o)), where v (f^) = Atr 6(0. — H), can be ex- 
pressed as a determinant 



Rn{{^a}) = detabKpfii^a^^b 



(2) 



In the Gaussian unitary ensemble (GUE) the kernel has 
the lirniting form Iuhn Kpf{fla,^b) = k{fla — ^b), 
where k {^la — ^b) = k {ila — ^b) — S {fla — ^b) [h — a) 
and, defining the average energy level spacing A, fc(51) = 
sin(7r57/A)/(7rf2/A). While the scaling variable A is sen- 
sitive to the specific properties of the system, the func- 
tion fc(f2) is universal provided Q, ^ £c, where is the 
time scale for ergodicity (equivalently, fi/A <C g, where 
g = £c/ A is the dimensionless conductance) |^ . 

Often one encounters applications where the interest 
is in the parametric dependence of the energy levels of 
a random Hamiltonian H' = H -\- V subject to an ex- 
tended non-uniform (but slowly varying) perturbation V 
such as that imposed by a magnetic field or by an exter- 
nal gate voltage on a quantum dot (QD). Surprisingly, 
such correlations exhibit the same degree of universal- 
ity: when expressed through the dimensionless param- 
eter — X"^ {{dei/ dXY) , where X parameterizes the 
strength of the external perturbation V (provided the 
distribution of H' is stationary with respect to X), the 



statistical properties of the entire random functions (x) 
become universal, dependent only on the fundamental 
symmetries of the system ^ . 

However, there exists a number of important applica- 
tions where the perturbation is either local, or rapidly 
fluctuating, or both. One example is presented by the 
orthogonality catastrophe in a disordered conductor ||] 
(of volume V) where, in the coordinate basis, V{y) = 
V V(5''(r— ro). Within random matrix theory (RMT), this 
is modeled by a rank r = 1 perturbation Vij = vNSijSn. 
The shape-distorting effect of a gate acting on a ballis- 
tic QD can be either global or local, depending on the 
relation between the spatial extent of the gate and the 
wavelength of the electron Xp. Similarly, attempts have 
been made to explore the pumping of electrons through 
a ballistic quantum dot under the action of two local 
gates ^ which effectively represent a rank r — 2 pertur- 
bation. The potential F(r) = u V[(5''(r - ri) - (5''(r - ra)] 
imposed by a bistable impurity in a disordered metallic 
host presents another application of a rank r — 2 per- 
turbation. A bistable dislocation segment provides an 
example of a perturbation of a still higher rank. The 
freedom to explore perturbations of high rank is easily 
afforded in experiments on microwave resonators with 
movable scatterers or in STM devices where the tip 
does not resolve the Fermi wavelength. 

The aim of this letter is to show that, when applied to 
a member of the unitary random matrix ensemble, per- 
turbations of arbitrary rank r can be incorporated into 
a general scheme which allows analytical expressions for 
the joint distribution function and the DoS correlation 
function to be determined explicitly. In particular, for 
a rank r matrix V — Nv whose non-zero eigenvalues 
form a set N{vk}'^^^, the (n, m)-point parametric cor- 
relation function between the initial and the final DoS 
i?„„({f]J,{f2'J) = {n:=iH^a)KL^i^'m can be 
expressed as the {n + m) X {n + m) determinant 



lim Rnrama}A^'b}) 



det 



kina-^b) V^},k{na-% 



(3) 
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Here k(na - n'^,) = fc(f^a ~ K') ~ S^ila - and 



T^v.n = det I 1 - v 



dn 



fc=i 



dr2 



(4) 



while the inverse of I?s,n has a convenient integral repre- 
sentation 

•' k=i 

The diagonal blocks of ^ reproduce the standard cor- 
relations (|^) (see also Ref. while, in the particular 
case of rank r = l, forn = TO=:l Eq. (^) recovers the 
result obtained in Ref. The corresponding joint dis- 
tribution function of the combined set of eigenvalues {ti\ 
of the matrix H and {e^} of H' is given by 



7'({6,},{6^})«7'({e.})||igldet,, 



4) 



(5) 



where 'P{{ei\) represents the distribution function of 
the eigenvalues of the matrix _ff for /3 = 2. 

To apply these results to physical systems, it is nec- 
essary to effect a generalization to account for cases 
where the Hamiltonian dynamics H is generated by a 
sum of regular and random components. This leads to 
non-vanishing off-diagonal terms in the ensemble aver- 
age propagator {G{€p)) = {{ep — In disordered 
metals, the corresponding non-universal terms are known 
as Friedel oscillations and reflect the short-ranged bal- 
listic nature of wave propagation via the Friedel func- 
tion /(r) = (AV/7r)Im(G(eF - iO;0,r)). In the pres- 
ence of an assembly of rank one perturbations V(r) = 
^j,WfcV^'^(r — Tfc), parametric spectral correlations are 
heavily influenced by the spatial arrangement of the 
wavefunction around the impurities. Surprisingly, the 
correlation function in this case retains the overall struc- 
ture of (^, and the effect of the Friedel oscillations is 
incorporated into the expression for Rnm by means of 
the generalization 



fv.n 



det 1 - fv 



d 

dn 



(6) 



where the elements of the r x r matrix / are prescribed 
by the Friedel function f^i = /(rj. — r;), and v is the 
diagonal matrix with entries VkSki- If the points and 
r; are separated by a distance in excess of the mean-free 
path, /(rfc — r;) i-^ Ski, and Eq. (Q) is recovered. 

Within RMT, Eq. applies to perturbations of ar- 
bitrary rank and strength. In systems where the Hamil- 
tonian dynamics includes a regular component (such as 
diffusive or ballistic conductors) Eqs. (Q) and (^) rep- 
resent the leading (zeroth) order term in an expansion 
in the inverse powers of the dimensionless conductance 
g. Whether the subleading terms in this expansion are 



small depends sensitively on the range and profile of the 
perturbation. If tr {fV)"^ > tr{VUV), where H is the 
diffusion (Perron-Frobenius) propagator in diffusive (bal- 
listic) systems, the results above apply without qualifi- 



cation 1 12, 1^. In the opposite limit, when rescaled by a 
single parameter which incorporates diffusive correlations 
at length scales in excess of the mean free path, para- 
metric correlations acquire a universal form Q , indepen- 
dent of the detailed properties of the system. The struc- 
ture of these correlations can be inferred from the RMT 
above: for an extended perturbation of rank r ~ 0{N), 
the requirement that V stays finite as — > oo implies 
that trV^ - ©(iV^A^). It follows that, for aU n > 3, 
tr{V/NA)" 0, and Vis,q exp[-{trv'^/2)d^/dn'^] 
reproducing identically the results of Ref. ||^ . 

In the standard Wigner-Dyson RMT, the determinan- 
tal structure of the many-point correlation functions 
makes it possible to obtain analytical results for proper- 
ties associated with individual energy levels, such as level 
spacing statistics and its various generalizations . Sim- 
ilarly, we anticipate that on the basis of the (n, m)-point 
parametric correlation function obtained in this Letter, 
an analogous theory of 'parametric level spacings' in the 
GUE can be developed. 

To illustrate this point, we consider the evolution of 
a single level under the influence of a rank r = 1 
perturbation. Specifically, we study the distribution 
P{d) = A{J2,S{n-e^)6[n + d- t[ {v)]) of distances 
^ ~ ^ between a given level Ci and its "descen- 
dant" e'i{v). This quantity is of independent interest in 
the context of Coulomb blockade (CB) peak spacings in 
disordered QDs. It is known that RMT-based statisti- 
cal theory of peak spacings shows poor agreement with 
experiment 1 14 . The disagreement is most likely at- 



tributable to the effects of strong Coulomb interaction 
which manifest themselves as "spectral scrambling" |l^ 

- readjustment of single-particle eigenstates to accommo- 
date the non-uniform spatial distribution of the charge of 
an added electron, and the parametric "gate effect" |l^ 

- distortion of the shape of the QD as the gate voltage 
is changed to sweep through a sequence of CB peaks. It 
would be interesting to study the latter effect in isola- 
tion, for example by comparing the two sequences of CB 
peaks corresponding to w = and v ^ Q. If the "gate 
effect" obeys RMT statistics, the distribution of differ- 
ences in peak positions should coincide with P{d). An 
r = 1 perturbation can be simulated, for example, by 
adjusting the voltage on a finger-shaped gate whose tip 
has a spatial extent smaller than Xp. 

A particular feature of the r = 1 perturbation is the 
identity < e- < e^+i where the levels are num- 
bered in increasing order (assuming w > 0). P{d) thus 
coincides with the distribution P{d), defined as the prob- 
ability that an interval of length d is bounded by a pair 
{ti,t'j(y)) at its ends, and does not contain inside it any 
levels from either of the two sets. The above identity 
ensures j = i. The standard method of computing level 
spacing distributions in RMT ||] affords a straightfor- 
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low at least two different routes: By implementing con- 
ventional methods of statistical field theory, Rnm can be 
expressed in the form of a supersymmetric field integral 
with an action of non- linear ti-model type |^ . Although 
such an approach has the advantage of allowing a sys- 
tematic development of non- universal corrections in both 
the disordered and chaotic environments, this method is 
largely tailored to the consideration of two-point corre- 
lation functions (i.e. n = m — 1). Thus, we will follow 
a complementary route employing a method based on 
orthogonal polynomials Q. A similar approach was em- 
ployed in Ref. to study the statistics of eigenvalues 
of matrices belonging to the GUE under the influence 
of a finite rank non-Hermitian perturbation. The latter 
ensemble serves as a model for describing resonances in 
quantum chaotic scattering JTot . 

The joint distribution of complex N x N Hermitian 
matrices H and H' is defined by 



FIG. 1: The P{d) = P{d) curves at several characteristic 
values of v for a rank r = 1 perturbation. 



V {H, H') = V{H)5{H' -H -V), 



(7) 



ward generalization to the parametric case allowing the 
distribution P{d) to be expressed as 

P{d) = det(l - K) [^11(0, 0)^22 (rf, d) - ^12(0, d)g2i(d, 0)] 
(cf. Eq. (5.4.31) of Ref. ||]), where the integral kernel 



K{e,e') 



kie-e') V^^lk{e-e') 
Ps,efc(e-e') fc(e-e') 



is defined on the segment [0, d] of the real line, and 
C/Q/3(e, e') is the inverse of 1 — K~^. In Fig. 1 we 
present a set of P{d) curves at several different val- 
ues of V. Note that in the limit t; ^ the distribu- 
tion Pv^o{d) = ie"*^/" coincides with the appropriately 
scaled level velocity distribution , while for v > 3 the 
curves P{d) develop a maximum at a finite value of d 
signifying weak repulsion between the old and the new 
levels. 

Let us turn now to the derivation of Eqs. (|3|), (Q), and 
(pf) . To obtain the correlation function Rnm one can fol- 



where V{H) represents the Gaussian distribution fmic- 
tion of 77, and the matrix (5-function is understood as a 
product of iV^ scalar (5-functions (one per each indepen- 
dent component of the complex Hermitian matrix H). 
Setting H = U7W and H' = C/'e' C/'^ where ? and ? 
represent the diagonal matrices of eigenvalues, the joint 
distribution function V {{ei} , {e'j}) is obtained from (|^) 
by integrating out the 'angular' degrees of freedom. 



/:^^{H)dii{u) 



\J(N) 



A^{H')dfi{U')V{H, H') 



U(W) 



where dji {U) represents the invariant Haar measure on 
the unitary group U(iV), and A{H) = A({eJ). To facili- 
tate the angular integration it is convenient to introduce 
a Lagrange multiplier in the form of the Hermitian ma- 
trix A = U\XUj^, where A represents the diagonal matrix 
of eigenvalues of A, after which the expression for the 
joint distribution function can be brought to the form 



^(fe},{ea)o^ / A'iH)dfi{U) A2(iJ')W)n/ A2(A)dAi(f/A) e'*'-[^(«+^-^')l-*'-^ 

JV{N) Jv(N) ,_i J -00 JV(N) 



IV{N) 



00 JlJ{N) 



The overall numerical constant can be restored by demanding the correct normalization. Integration over U and U' 
is performed by means of the Itzykson-Zuber integral pO| 

^„ (Tr\ (AUBU^) ^ de^ [exp (m^bj)] 
viNr^^ A(A)A(B) ' 
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where {a^}, and {hi} represent, respectively, the eigenvalues of the N x N Hermitian matrices A and B. As a result, 



^KX^iS) J\]{N) 1=1''-°° PP' 



N 



(8) 



where the sum over P runs over all permutations of the 
indices 1,...,N, and (—1)^ denotes the signature of the 
permutation. However, using the Itzykson-Zuber formula 
to undertake the remaining angular integration over U\ 
is impractical because of the N — r degenerate (zero) 
eigenvalues of the perturbation V. Instead, it is con- 
venient to express the rank r perturbation in the form 
V = N X)fc=i ''^k SLk (?) where a/c form a set of r com- 
plex mutually orthogonal A^-component vectors of unit 
length. Due to the invariance of the measure /i(t/), the 
integral in Eq. does not depend on afc, and there- 
fore its value does not change upon averaging over them. 
Although the measure dW of the integral over a^ must 
formally enforce the conditions of normalization and mu- 
tual orthogonality, in the large- A^ limit it can be replaced 
by the product of the Porter-Thomas distributions 
of their components: 



r N 



dw[{a}]=nn 



N daedal, _j^\^^^\2 



k=i 1=1 



Absorbing the rotations of A into the invariant measure 
dW, and integrating out the eigenvalues of A, one obtains 



V{{e,},m^n{e.}) 



A({q}) ^ 



N / 

I dW[{a}]n'5 4,(^)-ep(,)-A^^ 

i=l \ k=l 



Vk \aki\ 



(9) 



Then, introducing new variables Xki = V Naki , the inte- 
gral over Xki is identified as the integral representation 
of thus estabhshing Eq. (|). 

To obtain the correlation function i?„„i the joint distri- 
bution function (||) is multiplied by Yl"^i J2i^ =i ^i^a — 

^ia) rif^i Sjl=i ^(^6 ~ ^iJ ^^'^ then integrated over 
all eigenvalues {e,} and {e^}. To perform the in- 
tegrals over {ti} we utilize the following identity: 



nti2?~,',A({eJ) = A({e,}). Upon integrating by 
parts over all except the n fixed energies e^^ , and inte- 
grating over all , the correlation function can be recast 
as 



^ n \ m 



b=l 



(10) 



After some algebra, Eq. (Q) is straightforwardly estab- 
lished following the standard procedure of expanding the 
Vandermonde determinants in terms of Hermite polyno- 
mials and using their orthogonality properties. 

Finally, to establish the form of T) in the case when 

can be 



Friedel correlations are present, we note that Vij 

wp-»uz.fe=i -fe^i'Vi^^^lC^pV where 
the fc-th component of the ^-th eigenfunction of H . 



written as Y.iMY:k=i^k^t' i^k'^^Ul 



IS 

To 



the leading order in 1/g, the distribution of ■0^'^ is un- 
correlated with {ci}, and is given by Berry's conjecture 
p3| , In particular, in the unitary ensemble, 

dW [{V'}] oc J|dVi'^*dV'i'^e-S''^^'"(^"')..'^i' . 

kl 

Absorbing U into U\ and reversing the order of integra- 
tion over dW and U' , we arrive at Eqs. (^), (^. 

To summarize, we have demonstrated that the {n, m)- 
point parametric correlation function of the DoS in GUE 
possesses a determinantal structure for perturbations of 
arbitrary rank. Beyond GUE, a generalization of this 
determinantal structure is shown to accommodate spa- 
tial correlations induced by ballistic wave propagation at 
length scales shorter than the mean free path. 
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